Flexible Quasi-Three-Dimensional Terahertz Metamaterials 
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We characterize planar electric terahertz metamaterials fabricated on thin, flexible substrates 
using terahertz time-domain spectroscopy. Quasi-three dimensional metamaterials are formed by 
stacking multiple metamaterial layers. Transmission measurements reveal resonant band-stop be- 
havior that becomes stronger with an increasing number of layers. Extracted metamaterial dielectric 
functions are shown to be independent of the number of layers, validating the effective medium ap- 
proximation. Limitations of this approximation are discussed. 



Electromagnetic metamaterials enable numerous ex- 
otic responses not available in natural materials, such 
as negative refractive index perfect lensing and 
cloaking They typically consist of sub- wavelength 

metallic resonators fabricated on dielectric or semicon- 
ducting substrates and can be made to respond to either 
or both of the electric and magnetic fields of incident 
waves. By proper scaling, individual resonators and the 
composite metamaterial can resonantly respond to elec- 
tromagnetic waves of nearly any frequency, upto optical. 
Typically, the size of the resonators is roughly an order of 
magnitude smaller than the operating wavelength. Split- 
ring-resonators (SRRs) in combination with wire arrays, 
initially introduced by Pendry, have been widely used 
for metamaterial designs. These structures present a 
considerable fabrication challenge at higher frequencies 
(> 100 GHz) and therefore microwave techniques have 
dominated experimental metamaterial research. This is 
especially relevant for bulk, or three-dimensional (3D), 
metamaterials, where multi-layered structures become 
necessary. 

High frequency, multi-layer metamaterials can be fab- 
ricated by dicing and stacking single resonator arrays 
grown on thin substrates [4|,|5|. Or they can be mono- 
lithically grown on a single substrate using multi-layer 
processing [^, 0] ■ The latter process is particularly suit- 
able for optical metamaterials where the spacing between 
layers is only in the sub-micron or micron range. Either 
approach works for terahertz (0.1-10 THz) metamateri- 
als though the different spacing between successive lay- 
ers creates new challenges. An additional challenge, at 
all frequencies, is to create flexible or conformable meta- 
materials, which would be particularly useful in certain 
applications, such as cloaking or shielding. Fabrication 
techniques using rigid or fragile substrates, such as FR-4 
or silicon, are clearly unfit. However, some progress has 
been made here. Multi-layer metamaterials have been 
demonstrated in the far-infrared frequency regime by us- 
ing polyimide fillers between resonator layers . 

In the past few years, metamaterials research has 



sparked special interest in the THz community. Metama- 
terials are an optimistic approach to overcome the limi- 
tations of natural materials in the construction of novel 
functional THz devices, some of which have recently been 
demonstrated [8, 9., 10] ■ Planar SRR arrays fabricated 
on various semiconductor and insulator substrates have 
proven capable of responding to THz radiation electri- 
cally TlJ and/or magnetically [l^, [l3|- Other designs 
provide an app roach to negative permittivity Q or per- 
meability [12l | . 

In this work, we address the issues of flexibility, 3D fab- 
rication, and material parameter extraction by studying 
the transmission properties of THz metamaterials com- 
prised of metallic electric split-ring-resonators (eSRRs) 
on conformable polyimide (Kapton) substrates using ter- 
ahertz time domain spectroscopy (THz-TDS). These are 
studied as single layers and also as quasi-3D THz meta- 
materials by stacking up to four layers. In all cases, 
the measured transmission exhibits a clear minimum 
at 1.12 THz because of the excitation of the inductor- 
capacitor (LC) resonance. The transmission evolves into 
a full stop-band centered at the resonance frequency as 
the number of layers increases in the stacked quasi-bulk 
medium. Complex dielectric properties of the single and 
multiple layer samples are extracted and we show how 
they relate to the effective medium approximation. 

The eSRR structure employed in this work is shown in 
Fig. 1. The square resonators of length A = 40 /im are 
arranged in a square lattice of period P = 54 /im with 
line width = 3 /im, capacitor plate length I ~ 12 /im, 
and gap g — 3 /im. Using standard photolithography 
and e-beam evaporation, the eSRRs are fabricated in a 
square array on a commercially available polyimide film 
with a measured thickness of 84 ± 2 /xm. Polyimide 
has been widely used for photonic and electronic devices 
because of its high electrical and thermal stability [l3 |. 
Its flexibility, durability, relatively low refractive index 
{n = 2) and absorption (a = 20 cm^^) also make it 
favorable as a THz metamaterial substrate. Metalliza- 
tion consists of 10 nm of titanium followed by 200 nm of 
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FIG. 1: Diagram of split-ring-resonator unit cell. Black area 
represents the metallized resonator and gray area represents 
the unit cell. E and H shows the orientations of incident THz 
electric and magnetic fields. Bottom shows the schematic 
of the THz propagation through the metamaterial (region 3) 
sandwiched between the quartz plates (region 2). 



gold. During photolithography a regular silicon wafer is 
attached to the polyimide film to provide the mechan- 
ical support. The polyimide-based single layer meta- 
materials were then cut, visually aligned, and tightly 
stacked to form the quasi-3D media. All samples had 
a 10 mm x 10 mm patterned active area. 

THz transmission measurements were performed with 
wave incidence normal to the eSRR plane in a confocal, 
photoconductive antenna based THz-TDS system [l^. 
The system generates broadband pulsed THz radiation 
with a frequency independent beam waist of 3 mm diam- 
eter at the sample. Samples are placed in THz path with 
the gap bearing arm of the eSRR oriented parallel to the 
THz electric field, as shown in Fig. 1. The time-domain 
THz fields are recorded following passage through the 
samples and a dry air reference. All of the samples are 
sandwiched between two 1 mm thick quartz plates during 
measurements. Quartz has an index roughly equivalent 
to polyimide, which ensures that the front metamaterial 
layer has similar boundary conditions as subsequent lay- 
ers. The reference is measured with the quartz plates 
in contact with each other. The complex sample Es{lo) 
and reference Ej^{uj) spectra are calculated by numeri- 
cal Fourier transform. The effective material parameters 
of the sample can then be extracted from the measured 
transmission function t{Lu) = Es{lo) / Eji{lj) via [lij . 

exp {iknd{n - 1)) 
1 + f 237*32 exp(i2kodn) 

where ^23, ^32 and r23, are the frequency-dependent 
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FIG. 2: (a) Measured transmission amplitude of THz electric 
field through metamaterial samples with different number of 
layers, (b) The corresponding phase change. The measure- 
ments of 1-layer, 2-layer, 3-layer, and 4-layer metamaterials 
are represented by solid (black), dashed-dot (blue), dashed 
(red), and dotted (green) curves, respectively. 



complex transmission and reflection coefficients through 
the interfaces at the quartz-metamaterial boundaries; ko 
is the free-space wavenumber, d is the sample thickness, 
and h = n + IK is the sample complex refractive index. 
Time-windowing allows us to disregard multiple reflec- 
tions in the quartz plates. 

Figure 2(a) shows the measured THz amplitude fol- 
lowing transmission through the metamaterials. Mea- 
sured transmission spectra reveal clear resonances at 
1.12 THz due to the effective inductive-capacitive (LC) 
response of the eSRRs. Transmissions were also mea- 
sured through the multiple layer samples and then com- 
pared with the single layer case. The transmission spec- 
tra show band-stop behavior that increases in strength 
with an increasing number of layers. At four layers the 
effective medium becomes completely band-stop at the 
resonance frequency whereas the single layer transmis- 
sion minimum is 33%. The low frequency off-resonance 
transmission also decreases from 95% in a single layer 
metamaterial to 85% in a four-layer metamaterial. The 
overall linewidth of the LC resonance increases with num- 
ber of layers, though there is no significant shift in the 
resonance center frequency. This demonstrates that the 
resonance frequency is mainly determined by the SRR 
structural parameters. We also note the ripple features, 
most noticeable between the LC (1.12 THz) and dipo- 
lar (2.5 THz) resonances. Simulations initially indicate 
that these are due to the multiple reflections within or 
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FIG. 3: Extracted (a) real and (b) imaginary dielectric func- 
tions of the effective metamaterials. They are presented using 
the same patterns (colors) as in Fig. 2. The solid gray curve 
represents the dielectric function of the bare polyimide film. 



between sample layers. This is under investigation, but 
we presently restrict our discussion to the behavior of the 
LC resonance, our main interest. 

Figure 2(b) shows the phase changes incurred by the 
THz pulses during transmission through the metamateri- 
als. Both the substrate polyimide and the metallic eSRRs 
contribute to this phase change: the linear contribution 
comes from the substrate while the resonance structure 
comes from eSRRs. As expected, the slope of the lin- 
ear contribution increases as the wave encounters more 
substrate. 

The optical properties of the samples can be extracted 
from the measurements, but are critically dependent on 
the choice of some effective metamaterial thickness. Of- 
ten this effective thickness is defined by enforcing a cu- 
bic unit cell for each resonator 17[. In our case this 
would be 54 /im. Since our substrate is 84 /zm thick. 



this approach introduces complexity in the analysis and 
can even cause incorrect results [16]. Instead, we sim- 
ply assign each metamaterial layer the same thickness 
as our substrate, 84 fim. The complex dielectric func- 
tion is computed from the measured complex index by 
e{uj) = e,.(a;) + iei{uo) ~ v? . Figures. 3(a) and (b) show 
Cr and Ci, respectively, for the individual and stacked 
metamaterials, and the bare substrate. The bare poly- 
imide properties are fairly constant over our frequency 
range. The metamaterial properties are largely indepen- 
dent of the number of layers, in all the cases showing a 



fairly flat averaging about 3.8 with an approximately 
Lorentzian perturbation of about ±1 around resonance. 
The off-resonance metamaterial dielectric functions are 
very similar to those of the substrate with deviations 
that could be due to loss in the metal structures and 
tolerances in the substrate thickness. 



Figures 2 and 3 reveal how our metamaterials behave 
in light of the effective media approximation. Simply ex- 
pressed, €{lo) defines a material's electric response (e.g. 
phase shift) as a function of wave propagation distance. 
Thus, for homogenous media having fixed e(w), thicker 
samples obviously produce larger cumulative responses. 
Our multi-layer samples create larger responses than the 
single layer (see Fig. 2), but these are distributed over 
greater thicknesses (more layers), yielding an unchang- 
ing e{uj). This behavior is consistent with an effective 
medium approximation. Though not detailed here, this 
approximation fails at the individual layer. For any sin- 
gle layer the entire resonant response is produced by the 
very thin SRR array [l8| ; it does not accumulate in pas- 
sage through the substrate. For waves in transmission 
this fixed response is effectively distributed over the en- 
tire substrate thickness. Therefore, thicker substrates di- 
lute the response, and ultimately flatten the Lorentzian 
variation in e(u;). This behavior is obviously not char- 
acteristic of a truly homogeneous medium and marks an 
important limitation of the effective medium treatment 
when applied to planar metamaterials. 



In conclusion, we have demonstrated THz metama- 
terials based on eSRRs arrays fabricated on flexible, 
thin polyimide substrates and compared the transmission 
properties of the single-layer and quasi-3D bulk meta- 
materials using THz-TDS. Our results further verify a 
fabrication approach by which truly bulk 3D, durable, 
and conformable metamaterials can be realized; these 
are obviously important steps in realizing functional THz 
metamaterial devices such as prisms, lenses, waveguides, 
filters, and cloaks or shields. Finally, our results illus- 
trate interesting details about the proper application of 
the effective medium approximation as it applies to pla- 
nar metamaterials. This work was performed, in part, 
at the Center for Integrated Nanotechnologies, a U. S. 
Department of Energy, Office of Basic Energy Sciences 
nanoscale science research center jointly operated by Los 
Alamos and Sandia National Laboratories. Los Alamos 
National Laboratory, an affirmative action equal oppor- 
tunity employer, is operated by Los Alamos National Se- 
curity, LLC, for the National Nuclear Security Admin- 
istration of the U. S. Department of Energy under con- 
tract DE-AC52-06NA25396. We gratefully acknowledge 
the support of the U.S. Department of Energy through 
the LANL/LDRD Program for this work. 
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